Conditional sea surface statistics and their impact on geophysical sea surface parameters retrieved from SAR altimetry signals

In this study we present an analytical formulation of synthetic-aperture radar (SAR) altimetry signals including narrow banded nonlinear wave fields and conditional statistics between wave elevation displacements, horizontal wave slopes and vertical wave particle velocities. Considering the wave elevation displacements coskewness with respect to horizontal slopes leads to an analytical formulation of the electromagnetic bias within a SAR-mode altimeter stack. This formulation can be either parametrized by the significant wave height (SWH) and mean wave steepness, or in terms of the variance of vertical wave velocities. The effect of conditional vertical wave particle velocity variances with respect to the observed horizontal wave slopes close to nadir incidence angles leads to an effective reduction of the azimuth blurring of SAR-mode stacks. We present here a formulation of this effect by examining JONSWAP ocean wave spectra. In most cases this effect reduces the azimuth blurring by 10% to 30%. Additionally we investigate the effect of a nonlinear wave elevation displacement probability density function (PDF) on estimated geophysical parameters. We were able to show that including an elevation displacement skewness of 0.13 improves significantly the SWH consistency between altimetry and ECMWF Reanalysis v5 ERA5 results.All of these effects are validated with respect to ERA5 model data in the North East Atlantic region and in situ data located in the German Bight and Baltic Sea.The developed model can be used in both SAR and conventional altimetry retrackers.     

详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

复杂动态环境下无人飞行器动态避障近似最优轨迹规划

针对复杂动态环境下无人飞行器的动态障碍规避问题，基于合理假设建立了无人飞行器和动态障碍的运动学模型，并综合考虑无人飞行器飞行过程中的终端约束、控制输入约束、安全避障约束等，以能量最少为性能指标构建动态避障问题数学描述。之后，针对终端约束和控制输入约束，依据优化模型预测静态规划算法(OMPSP)生成初始轨迹；针对动态避障问题的不等式约束，引入松弛变量并结合滑模变结构控制方法设计松弛变量动力学，实现对一个、多个或同时多个动态障碍的安全规避；最后，依据有限时间微分动态规划(RHDDP)算法进行轨迹优化，获得满足上述各种约束并能规避动态障碍的近似最优轨迹。     

基于双滤波器的高精度目标跟踪估计方法

针对连续推力的合作航天器，采用双重无迹卡尔曼滤波(DUKF)算法估计其状态和加速度。通过状态滤波器和参数滤波器的配合，提升滤波精度，完成运动状态和参数的估计，从而实现合作目标的运动轨迹跟踪。与合作航天器相比，非合作航天器存在大小未知、发生时刻未知的机动，无法获得加速度，且信息获取和运动状态的估计难度大。针对非合作航天器，基于简化的相对运动方程，结合天基平台获得目标的观测信息，采用两个扩展卡尔曼滤波(EKF)及基于半正焦弦的机动检测策略实现多未知脉冲机动的运动状态的估计。仿真结果表明：相比于无迹卡尔曼滤波（UKF），DUKF在对合作航天器的状态和加速度估计方面具有更快的收敛速度和更高的滤波精度；对于存在未知机动的非合作航天器，通过对比验证机动检测策略与滤波器切换策略相结合的方法的有效性，该方法能够检测到多次机动并且减少误判。     

星载激光通信终端在轨机动对温度影响

对激光通信终端在轨瞬态温度变化开展了仿真，以期研究在轨机动的影响、热分析和热试验时机动模式的简化模拟。通过合理地分析与简化，建立终端的轨道热分析模型，准确模拟在轨机动，根据典型机动模式、外热流和涂层退化等因素设计了计算工况，得出了终端在轨运行过程中的温度场随时间和姿态变化规律。结果表明：不同机动模式下的温度变化存在差异，最大可达23.0℃，采用固定姿态或一维转动的简化热分析或热试验不能准确模拟实际飞行温度，甚至不能部分替代二维转动热分析或热试验；光学天线和反射镜的温度控制是制约终端工作的瓶颈，为终端设计合适的避光机动策略，可大幅度提高温度稳定度和均匀度。研究结果可以为光机电设备在轨机动策略的设计和改进提供参考。     

轨控大干扰下的系统角动量管理

摘要： 针对利用角动量交换装置系统实现轨控期间姿态维持的情况，研究当轨控推力器存在大力矩扰动时系统角动量管理方法.结合对地运行航天器轨道运动特性，分析轨控干扰力矩的积累角动量变化规律；考虑角动量管理装置系统的角动量存贮容量，基于角动量积累规律提出了一种结合偏置角动量建立与对称分布轨控位置选择的系统角动量管理方法，实现了对称轨控位置积累角动量互相抵消的自平衡效果，极大程度地提高轨控效率且避免系统角动量饱和现象.所提出的方法适用于喷气欠驱动控制系统，方法的有效性通过了数学仿真验证及在轨型号应用.     

拦截高速目标的比例与反比例导引捕获区分析

针对拦截高速目标的作战特点，分析了比例导引(PN)与反比例导引(RPN)的捕获区。首先，通过分析拦截弹与目标的相对运动关系，推导得到了顺轨和逆轨的零控拦截条件，此条件由目标和拦截弹的速度前置角以及二者速度比确定；其次，以拦截弹和目标速度前置角为坐标系，推导得到了PN以及RPN捕获区以及各自导航比设置范围。PN的捕获区由逆轨零控拦截条件以及与其相切且斜率为1/(N-1)的两条直线构成，RPN的捕获区由顺轨零控拦截条件以及与其相切且斜率为1/(-N-1)的两条直线构成；然后，利用函数对称性将PN与RPN捕获区转换到同一坐标区间，得到了相同条件下RPN捕获区要大于PN捕获区的结论；最后，开展了四种情形下的仿真，验证了本文捕获区分析的合理性及有效性。     

不同密度弹丸对水冰的超高速撞击成坑实验

为了研究冰冻天体表面撞击坑的形成与演化，开展了水冰的超高速撞击成坑实验。使用二级轻气炮发射1.0 mm直径的球形弹丸，以3 km/s、5 km/s和7 km/s速度对圆柱状冰块进行撞击。弹丸材料包括聚碳酸酯和不锈钢两种，冰块温度为253 K。实验观察到了不同弹丸和不同速度条件下，冰块中撞击坑的形貌特征。对撞击坑直径、深度和剖面形状进行了测量，并与文献中铝弹丸对水冰的撞击坑进行了比较分析。获得了水冰撞击坑特征随撞击参数的变化规律，结果表明：撞击坑直径和深度的主导机制不同，坑深主要由弹丸侵彻作用形成，而坑径主要由冰块的剥落所致；坑深比坑径具有更强的对于弹丸密度的依赖性，高密度弹丸撞击坑直径具有比低密度弹丸更强的对于撞击速度的依赖性；撞击坑体积与撞击能量成正比，高密度弹丸形成的撞击坑直径表现出“能量缩比”行为，而低密度弹丸形成的撞击坑直径表现出“动量缩比”行为。     

椭圆轨道卫星对地全球覆盖电推进控制

针对椭圆轨道卫星近/远地点的星下点对全球或特定纬度区域的访问问题，提出一种连续小推力下的对地覆盖控制策略。首先，推导了自然摄动对卫星拱线变化的影响，并探讨了进行小推力覆盖控制的必要性。然后，针对燃料消耗的优化问题，将控制方程展开成含傅里叶级数的形式，用以获得便于星上计算的解析形式的次优解，同时探讨了截取阶数与优化程度的关系。在进行拱线控制的同时，通过合理设置约束，对椭圆轨道的近地点高度进行保护，确保卫星安全运行。仿真结果表明，提出的方法能够以适当的燃料消耗代价实现椭圆轨道的近/远地点的全球覆盖控制或特定纬度区域的反复推扫，且控制力在可接受的范围内。     

Substitution of satellite-based land surface temperature defective data using GSP method

Land surface temperature (LST) as an important environmental variable provides valuable information for earth environmental system modelling. Currently, LST is obtained through satellite thermal sensors at various spatial and temporal resolutions. Although spatially continuous satellite-based LST measurements are intended to overcome the shortcomings of sparse ground-based LST measurements, LST images often contain anomalous values due to the existence of clouds or sensor malfunctioning. The problem becomes more serious where the users deal with high spatial resolution characterized by low temporal resolution. This study examines the capability of a newly developed graph signal processing (GSP) method using two-dimensional single-date thermal data. For this purpose, four Landsat/TIRS datasets are analyzed. The data of five elliptical regions on thermal images are eliminated and then reconstructed through the GSP method and using the LST values of the enclosing rectangles containing the ellipsoids. The results indicate that the temperature variation determined by the GSP method generally conforms to the original image LST values. According to a correlation test conducted on the original image LST and those obtained through the GSP method, the values vary from 58% to 95%, which is an above-the-average rate (RMSE from 0.69 to 2.27). The statistical analysis of the original image LST in both the elliptical regions and the enclosing rectangles containing the ellipsoids indicates that an increase in the variance of LST data causes an increased error in the calculation of temperature by the GSP method, and vice versa. The results of the analysis of variance (ANOVA) and Duncan test indicated that an increase in the number of the non-zero spectral bins would result in increased RMSE values for all the dates and the regions. Moreover, the model errors were significant at the 0.05 level across all the image date and five elliptical study regions. Based on the results, the use of this method is recommended for the reconstruction of LST missing values, where dissimilarity of atmospheric conditions limits the use of other methods that depend on the time series data of various dates and a great deal of data calculation.